Superconducting properties of MgB2 from first principles 
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Solid MgB2 has rather interesting and technologically important properties, such as a very high 
superconducting transition temperature. Focusing on this compound, we report the first non-trivial 
application of a novel density-functional-type theory for superconductors, recently proposed by 
the authors. Without invoking any adjustable parameters, we obtain the transition temperature, 
the gaps, and the specific heat of MgB2 in very good agreement with experiment. Moreover, our 
calculations show how the Coulomb interaction acts differently on a and tt states, thereby stabilizing 
the observed superconducting phase. 
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Understanding and predicting superconducting prop- 
erties of real materials is of both fundamental and tech- 
nological importance. While there are interesting classes 
of materials, such as the high Tc Cu oxides, where the 
superconducting mechanism is still under debate, the re- 
cent discovery of phonon-mediated superconductivity at 
39.5 K in MgB2m - but also in other materials - keeps 
the interest in the phonon-driven mechanism alive. MgB2 
has rather peculiar properties, such as the presence of 
two superconducting gaps at the Fermi level. While two- 
band or, more generally, multi-band superconductivity 
has long been known to favor a high critical tem- 
perature, there remains the challenging question: Gould 
anyone have predicted quantitatively the peculiar super- 
conducting phase in MgB2, including its high Tc value, by 
means of "physically unbiased" , first-principles calcula- 
tions? 

As a matter of fact, despite the theoretical and tech- 
nological interest involved, the capability to predict from 
first principles material-specific properties, such as the 
critical temperature and the gap, has been out of reach 
so far. This is because (conventional) superconductiv- 
ity appears as the result of a subtle competition between 
two opposite effects, the phonon mediated attraction (de- 
noted "e-ph" in the following) and the direct Coulomb 
repulsion (denoted "e-e") between the electrons. His- 
torically, after the microscopic identification of the su- 
perconducting order parameter by Bardeen, Gooper and 
Schrieffer (BGS) [3, the first theoretical framework aim- 
ing at the description of real materials was put forth 
by Eliashberg In this theory the e-ph interaction 
is perfectly accounted for, however the effects of the e- 
e Goulomb repulsion are condensed in a single parame- 
ter, II* , which is difficult to calculate from first principles 
and which, in most practical applications, is treated as 



an adjustable parameter, usually chosen as to reproduce 
the experimental Tc. In this sense, Eliashberg theory, in 
spite of its tremendous success, has to be considered a 
semi-phenomenological theory. 

Looking at normal-state properties, density functional 
theory (DFT) 5] has enjoyed increasing popularity as a 
reliable and relatively inexpensive tool to describe real 
materials. In 1988 the basic concepts of DFT were gen- 
eralized to the superconducting phase by Oliveira, Gross 
and Kohn (OGK) via including the superconducting 
order parameter as an additional "density" in the for- 
malism. The theory of OGK, however, was restricted to 
weak e-ph coupling. Along the lines of a recently pre- 
sented multi-component DFT of electrons and ions, 
the theory of OGK was successfully generalized to the 
strong-couphng case 8]. In this Letter we apply the 
strong-coupling DFT to the challenging case of MgB2. 

The central equation of the DFT for superconductors 
is a generalized gap equation of the form 

1 tanh (^|S„'k') 

^nk — ^-2„k A„k 77/, ^nk,rt'k' ^n'k'j 

^ ■TT, I^n'k' 
n'k' 

(1) 

where n and k are respectively the electronic band index 
and the wave vector inside the Brillouin zone. (3 is the in- 
verse temperature and £'„k = \J (cnk — m)^ + l^nkl^ are 
the excitation energies of the superconductor, defined in 
terms of the gap function A„k, the Kohn-Sham eigenen- 
ergies of the normal state e„k, and the chemical potential 
fi. The kernel, /C, appearing in the gap equation consists 
of two contributions /C = IC^^^^ + /C°~°, representing the 
effects of the e-ph and of the e-e interactions, respec- 
tively. The first of these terms involves the e-ph coupling 
matrix, while the second contains the matrix elements 
of the screened Goulomb interaction. Eq. |^ has the 
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same structure as the BCS gap equation, with the kernel 
K, replacing the model interaction of BCS theory. This 
similarity allows us to interpret the kernel as an effec- 
tive interaction responsible for the binding of the Cooper 
pairs. On the other hand, Z plays a similar role as the 
renormalization term in the Eliashberg equations. We 
emphasize that Eq. is not a mean-field equation (like 
in BCS theory), since it contains correlation effects. Fur- 
thermore, it has the form of a static equation - i.e., it does 
not depend explicitly on the frequency - and therefore 
has a simpler structure than the Eliashberg equations. 
However, this certainly does not imply that retardation 
effects are absent from the theory: as a matter of fact, 
an Eliashberg-type spectral function a'^F{Q) enters the 
calculation of Z and A;;°"P^. 

Coming back to MgB2, its Fermi surface has several 
sheets with different orbital character (see e.g. Ref. 0). 
In particular, the tubular structures with a character are 
very strongly coupled to the E2g phonon mode, corre- 
sponding to a B-B bond-stretching in the boron planes 
MgB2 also has three-dimensional tt bands, that 
give rise to a complicated Fermi surface. Without holes 
in the a bands, the compound would not be supercon- 
ducting, like AIB2. The tt bands are coupled much less 
efficiently to phonons, but are nevertheless crucial to 
superconductivity. A remarkable feature of this com- 
pound is the presence of two gaps on the u and tt bands, 
as clearly demonstrated b y se veral different experiments 
d in 113, III III El 113 El- On the theoretical 
side, this system has been treated within the k-resolved 
Ehashberg theory 0,|2^, using a two-band scheme [2lj |. 
with four e-ph spectral functions to represent the dis- 
tinct couplings. Correspondingly, the anisotrop y o f the 
Coulomb interaction has also been investigated^ I^^l , 
with ^* treated as a matrix. 

In our density functional calculations we used the four, 
band resolved, Eliashberg functions, {a^Fn^n'i^y^n^n' = 
cr, tt), previously employed within a two-band Eliashberg 
scheme by Golubov et al. IH]. Our procedure keeps the 
fundamental distinction between a and tt gaps, analo- 
gously to the Eliashberg calculations reported to date. 
We recall here that a fully consistent calculation should 
not use the a^F„_„/ (fi) functions, but rather the nk, n'k'- 
resolved coupHngs: By using the a^-F„_„' (17), the e-ph 
interaction is averaged over k and k' at the FS, which 
may have a small, but non-negligible effect. We define 
the bands crossing the FS to be of a character if they 
are contained in a cylinder of basis radius 1/4 of a recip- 
rocal lattice vector, and of tt character otherwise. Away 
from the Fermi surface this distinction is meaningless but 
harmless, as the phonon terms die off quickly and the 
Coulomb term is independent of this distinction. As we 
described our computational strategy elsewhere 0, we 
shall not report the details here. We just mention that 
extreme care needs to be taken with the sampling of the 
region in k-space around the FS. We accomplish this task 
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FIG. 1: Calculated superconducting gap of MgB2 as a func- 
tion of energy (T = K). 

by using around 65000 independent k points, chosen ac- 
cording to a Metropolis algorithm. This places ~ 8000 
points within « tiLOE^g of the FS. We estimate the overall 
numerical stability of our results to be roughly 5%. 

As for the e-e interaction, we calculated the matrix 
elements of the screened Coulomb potential with respect 
to the Bloch functions of MgB2, for the whole energy 
range of relevant valence and conduction states, using the 
full-potential linearized augmented plane wave method. 
In our previous work the Coulomb interaction was 
screened using a simple Thomas-Fermi model. However, 
the different nature of a and tt bands in MgB2, and in 
particular the highly localized character of the former, 
strongly calls for the use of a non-diagonal screening, 
including local field effects. To avoid the cumbersome 
calculation of the dielectric matrix of MgB2, we used a 
dielectric matrix obtained with the method of Hybertsen 
and Louie 25] . In particular, we substituted the model 
dielectric function in Eq. (7) of Ref. IH by a Thomas- 
Fermi model, computed at the local (valence) density of 
MgB2. The diagonal part of our model cornpares rather 
well with the calculations of Zhukov et al. [2y|. 

In Fig. n we plot the energy gap as a function of (pos- 
itive) energy distance from the Fermi energy (the gap 
function exhibits particle-hole symmetry to a good ex- 
tent). We can see that the two gaps of the material, A^- 
and Att, appear naturally from our calculations. The tr 
gap is defined only up to the energy of the top of the 
a band. Both Ao- and A^r are anisotropic. This results 
from the anisotropy of the Coulomb potential matrix ele- 
ments - roughly 0.4 meV, « 6% of Acr at the FS and gets 
much larger at high energy, where there are many bands 
with different characters. The averages of A^ and Att 
at the Fermi level (6.8 meV and 2.45 meV respectively) 
are in excellent agreement with experiment. Both gaps 
change sign, which is a necessary condition to find super- 
conductivity in the presence of the repulsive Coulomb 
interaction. In fact, our gap equation does not converge 
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FIG. 2: Superconducting gaps at the FS and specific heat of 
MgB2 . Panel (a): Comparison between theoretical and ex- 
perimental gap at the FS plotted as a function of temperature. 
The calculated gaps and Tc (34.1 K) are obtained without the 
use of any adjustable parameter. Panel (b): Experimental 
and calculated electronic specific heat, as a function of T/Tc. 



to a superconducting solution, unless all electronic states 
in a large energy window are included. 

In Fig. |2] (panel (a)), the superconducting gaps are 
plotted versus temperature, together with a few recent 
experimental results. The agreement is striking: the val- 
ues of Tc (34.1 K) and of and A^^ at T = K are very 
close to the experimental data. Moreover, the temper- 
ature behavior of both gaps, along with their strongly 
non-BCS behavior, are very well reproduced. We believe 
that such an agreement for a highly non-trivial supercon- 
ductor such as MgB2, without any adjustable parameter, 
is unprecedented in the field of superconductivity. Ob- 
viously, unlike calculations performed using Eliashberg 
theory, we do not reproduce exactly the experimental 
critical temperature, as our calculations are not fitted to 
match any experimental quantity. 

We also calculated the Kohn-Sham entropy as a func- 
tion of temperature and, from its temperature derivative, 
the specific heat. In order to compare our results with 
experiments US 113, Ei, we plot in Fig. El (panel (b)) the 
reduced specific heat versus temperature, normalized to 
Tc (using the corresponding experimental and calculated 
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FIG. 3: Electronic repulsion contributions to the gap equation 
(see Eq.0, at an arbitrary k, averaged over k' and bands on 
isoenergy surfaces, but keeping the a, tt distinction. 



Tc values). Both the shape of the curve as well as the 
discontinuity at Tc are almost perfectly reproduced. We 
recall that the low temperature shoulder comes from the 
presence of the smaller tt gap and that our Aa-fA^^ is 
slightly different from the experimental ratio. 

While the good agreement with experiment underlines 
the predictive power of our method, it is only one part 
of our investigation. Another important aspect is to 
gain further insight into the peculiar superconductivity of 
MgB2. To this end, we performed a calculation using an 
average functional form for the Coulomb interaction, not 
distinguishing between a and tt bands. This functional, 
described in detail in Ref. |^, corresponds to a semiclas- 
sical treatment valid in the limit of slowly varying den- 
sities, and actually leads to a good agreement with the 
full matrix element calculation for s,p- metals, and works 
reasonably also in d-metals as Nb. We obtain Tc =52 K, 
with the a and tt gaps being 9.8 meV and 1.9 meV re- 
spectively. This test shows that the repulsion among a 
states, stronger than within tt and between a and tt, is 
crucial in achieving good agreement with experiment. We 
also see clearly that the more delocalized tt electrons are 
reasonably well described by an average formula derived 
from free-electron concepts, while the repulsion among 
a electrons needs the real matrix element calculations. 
Reversely, if we average all the e-ph spectral functions, 
we obtain a single gap A=3.7 meV and the much lower 
Tc ~ 20 K, in agreement with the analysis of Ref. 
and with a similar test carried out by Choi et al. |l9l | . 

Of course, in the search for novel superconductors with 
higher transition temperatures, it would be desirable to 
keep the extremely strong e-ph coupling for a states, 
while reducing the corresponding Coulomb interaction to 
that of more delocalized tt states. Unfortunately, the two 
features are linked together, as both the strong e-ph cou- 
pling and the strong Coulomb repulsion derive from the 
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covalent nature of the corresponding electronic states. 

To push this analysis further, we plot, in the up- 
per panel of Fig. 13 the energy dependence of the 
Coulomb contribution to Eq.(Q at T w OK, namely 

^™k,n'k'^n'k'/(2i^n'k') = ^^'k,„'k' Xn'k' agaiust ^n'k' = 

e„'k' — Ep, for a few k-points arbitrarily chosen on the 
Fermi surface. As Xn'k' goes to 0.5 on the Fermi surface. 
Fig. 131 shows the larger magnitude of the intraband a-a 
and TT-TT, relative to the interband a-n matrix elements 
(see also Ref. Esj) . The different behaviour of the tt 
and TT— > a terms results from the x factor, as the matrix 
elements themselves are symmetric. Obviously, the a-a 
repulsive matrix elements are the strongest. The scatter- 
ing of data for a given energy comes from the different 
orbital character of wavefunctions at n'k'. The energy 
dependence of the quantities plotted in Fig. Incomes al- 
most entirely from Xn'k', as the matrix elements them- 
selves have a marginal energy dependence. In order to 
show how the reduction of Coulomb repulsion takes place, 
we emphasize in the lower panel of Fig. |31the region cor- 
responding to low matrix elements (on a linear scale). 
We see how, when A becomes negative, the Coulomb in- 
teraction actually gives a constructive contribution due 
to the minus sign in Eq. Although the negative val- 
ues are much smaller (by almost 3 orders of magnitude) 
than the positive ones, the corresponding energy range is 
much larger, resulting into the well known, substantial re- 
duction of the effective Coulomb contribution. The most 
important energy region is located in between 0.50 to 3 
eV below Ep, in particular due to the strong intraband 
a-a matrix elements. The interband contribution from tt 
bands (violet in Fig. (S)), on the other hand, is consider- 
ably smaller, which is obviously the case also for the a 
contribution to (green in Fig. ^ . Summing up over 
n'k', the negative contribution to A^- coming from the 
(positive gap) region of the a Fermi surface overcomes 
by a factor of w 7 the contribution coming from the cor- 
responding TT region. 

Finally, it is also interesting to investigate the impor- 
tance of local field (LF) effects on the superconducting 
properties of MgB2. It turns out that using a diago- 
nal Thomas-Fermi (TF) screening (which completely ne- 
glects LF's) the a-a matrix elements increase by roughly 
15% (red in Fig. OJ, while the a-ir and tt-tt terms re- 
main basically unchanged. As a consequence, neglecting 
LF effects leads to a marginal (w 4%) decrease of A^r, 
but decreases significantly A^. (by about 14%). The dif- 
ferent behavior of and A^ can be understood quite 
naturally: LF corrections imply a better screened inter- 
action among electrons when they are located in a high 
density region inside the unit cell. This is precisely the 
case of the a bands. On the other hand, the tt bands 
are more delocalized - the electrons reside in the inter- 
stitial region - and are therefore reasonably described by 
diagonal screening. 

In this communication we presented the first non- 



trivial application of a recently developed ab-initio the- 
ory of superconductivity. In particular, we obtained for 
MgB2 the value of Tc, the two gaps, as well as the spe- 
cific heat as a function of temperature in very good 
agreement with experiment. We stress the predictive 
power of the approach presented: being, by its very na- 
ture, a fully ab-initio approach, it does not require semi- 
phenomenological parameters, such as ji* . Nevertheless, 
it is able to reproduce with good accuracy superconduct- 
ing properties, up to now out of reach of first-principles 
calculations. Furthermore, our calculations allow for a 
detailed analysis of the contribution of the Coulomb re- 
pulsion to the superconducting gap, opening the way to 
tailoring the electronic properties of real materials in or- 
der to optimize superconducting features. 
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